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ABSTRACT 

The properties of the interstellar medium (ISM) surrounding a planetary system can impact plane¬ 
tary climate through a number of mechanisms, including changing the size of the astrosphere (one of 
the major shields for cosmic rays) as well as direct deposition of material into planetary atmospheres. 
In order to constrain the ambient ISM conditions for exoplanetary systems, we present observations 
of interstellar Na i and K i absorption towards seventeen early-type stars in the Kepler prime mission 
field of view. We identify 39 Na i and 8 K i velocity components, and attribute these to eleven ISM 
clouds. Six of these are detected towards more than one star, and for these clouds we put limits on 
the cloud properties, including distance and hydrogen number density. We identify one cloud with 
significant (> 1.5 cm“^) hydrogen number density located within the nominal ^100 pc boundary of 
the Local Bubble. We identify systems with confirmed planets within the Kepler field of view that 
could lie within these ISM clouds, and estimate upper limits on the astrosphere sizes of these sys¬ 
tems under the assumption that they do lie within these clouds. Under this condition, the Kepler-20, 
42, and 445 multiplanet systems could have compressed astrospheres much smaller than the present- 
day heliosphere. Among the known habitable zone planet hosts, Kepler-186 could have an astrosphere 
somewhat smaller than the heliosphere, while Kepler-437 and KOI-4427 could have astrospheres much 
larger than the heliosphere. The thick disk star Kepler-444 may have an astrosphere just a few AU 
in radius. 

Subject headings: ISM: clouds — ISM: lines and bands — ISM: structure — local interstellar matter 


1. INTRODUCTION 

The interstellar medium (ISM) is often studied by ob¬ 
serving absorption lines in stellar spectra, which have 
been imprinted by intervening interstellar matter (e.g., 
lHartmaimlll904H . For ground-based optical spectra this 
can most easily be accomplished by observing early-type 
stars, which have spectra largely free of stellar absorp¬ 
tion lines; any stellar lines that are present are typically 
significantly rotationally broadened. Interstellar matter 
imprints a signature of narrow absorption lines on the 
stellar spectrum. High resolution spectra allow the sep¬ 
aration of multiple velocity components (corresponding 
to different ISM clouds), typically separated by a few km 
s“^. The primary ions accessible in the optical are Na i, 
K I, and Ca ii, while many other species are ob servable in 
the ultraviolet le.g.. iR.edfield fc Linskvl[^04ajl . The dis¬ 
tance to the star provides an upper limit on the distance 
to the interstellar cloud, while the most distant stars in 
the same region of the sky that show no absorption pro¬ 
vide a lower limit to t his distance fe.g.. iLallement et al.l 
l2003l[Peek et al.l[2011h . This technique has been used to 
map the three-dimensional structure of the ISM on both 
large and small scales. For instance, ILallement et al.l 
(|2003f) mapped th e ISM in all directions out to distances 
of ^ 250 pc, while iRedfield &: Linskvl (1200811 mapped the 
structure and kinematics of the local ISM (LISM), the 
collection of clouds in the immediate vicinity of the solar 
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system. Most previous work using these techniques cov¬ 
ered much or all of th e sky (e.g., iWeltv fc Hobbs! 120011 : 
ILallement et al.l l2003|l . Targeted surveys of a few se¬ 
lected fi elds have been conducte d, such as towards the 
Hyades (|Redfield fc Linsk^l2001ll an d in the direction of 
the s olar system’s past trajectory (jWvman fc Redfi^ 
1201311 . In this paper we present the first targeted survey 
of the ISM in a direction of exoplanetary interest: the 
Kepler field of view. 

In order to provide motivation and background for this 
work, in gru we introduce the interactions between the 
heliosphere (and by extension astrospheres) and the ISM, 
while in lH.2l we review the effects of variable astrosphere 
size upon any planets residing within. In H.dl we discuss 
how our observables in terms of the ISM and exoplan¬ 
etary properties can be used to constrain astrospheric 
properties, and provide a brief introduction to the Ke¬ 
pler mission. 

1.1. The Heliosphere and the ISM 

In our own solar system, the solar wind carves out 
a cavity in the surrounding ISM known as the helio¬ 
sphere. As the solar wind moves outward, it eventually 
encounters the ISM, resulting in a series of shocks and 
boundaries. The innermost of these is the termination 
shock, where the solar wind decelerates from supersonic 
to subsonic velocities. Both Voyager 1 and Voyager 2 
have crossed this shock, at distances of 94 and 84 AU, 
respectively (jStone et ai.ll2005l 1200^ . The correspond¬ 
ing outer boundary where the ISM becomes subsonic, 
the bow shock, is expected to lie at a distance of ^ 250 
AU (|Izmodenov et all 1200511 : however, recent work has 
suggested that the solar system is moving through the 
LISM at an ins ufficient velocity to ind uce the formation 
of a bow shock (|McComas et al.ll20I2ll . In between these 
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two shocks lies the heliopause, the contact interface be¬ 
tween the solar wind and the LISM. Voy ager 1 crossed 
the h eliopause at a distance of 121 AU (|Gurnett et alJ 
1201,111 ■ while Voyager 2 has yet to cross. The turbulent 
region between the termination shock and the heliopause 
is known as the heliosheath. Other stars also have stellar 
winds, and so will possess analogues to the heliosphere, 
known as astrospheres, with analogous structures (i.e., 
termination shock, astropause, bow shock, astrosheath). 

Planetary surfaces can be protected from cosmic rays 
by three major shields-the interplanetary magnetic field 
carried by the solar (stellar) wind within the heliosphere 
(astrosphere), the planetary magnetic f ield, and the plan¬ 
etary atmosphere (iScherer et al.ll200CT . Because the so¬ 
lar (stellar) wind carries its embedded magnetic field as 
far as the heliopause (astropause), when we refer to the 
“heliosphere size” or “astrosphere size” for the remainder 
of this article, we are referring specifically to the distance 
from the star to the heliopause or astropause. We note 
that the heliosphere and astrospheres are not spherical, 
but are rather swept back into a “cometary” shape by 
the star’s motion through the ISM; the sizes that we 
quote are the distances to the astropause in the upwind 
direction. The importance of the heliopause is illustrated 
by the fact that the cosmic ray flux measured by Voy¬ 
ager 1 increased by ~ 9% in less than a d ay when it 
crossed the heliopause (|Krimigis et al.l[2M3l) . The mod¬ 
ulation of cosmic rays within the heliosphere is caused 
primarily by cosmic rays scattering through magnetic 
fields in the turbu lent solar wind and the heliosheath 
(|Muller et al.ll2006ll . Thus, the efficacy of the heliosphere 
as a cosmic ray shield will be affected by both the size 
of the heliosphere and the degree of turbulence. The 
turbulence is caused in large part by th e scattering of 
pickup ions within the heliosp here (e.g., Ilsenberg et al.l 
I200.H iGamavunov et al.l 120121 1. which is in turn related 
to the neutral density in the surrounding ISM. 

The size of the astrosphere surrounding a star is de¬ 
termined by a momentum balance between the outgo¬ 
ing stellar wind and the streaming ISM. Higher ISM 
momentum-i.e., higher ISM density or higher yelocity 
of the star relative to the ISM-will reduce the size of 
the astrosphere. If the size of the astrosphere is reduced, 
the distance over which cosmic rays cross magnetic field 
lines is reduced, and thus the shield is less effective; the 
cosmic ray flux at any pla nets will increase. Quantita¬ 
tively, iMiiller et al.l (j2006[ ) found that during a passage 
through a dense interstellar cloud where the heliopause 
lies at a distance of 26 AU from the sun, the intensity 
of cosmic rays with energies between 300 MeV and 1 
GeV at the Earth will increase by a factor of 1.4 to 
7.6 over the present flux, depending upon assumptions 
about the power spectrum of the solar wind turbulence, 
corresponding to the relative importance of cosmic ray 
modulation in the solar wind versus in the heliosheath. 

The ISM in the immediate vicinity of the solar system, 
the LI SM, has been extensively studied; see lFrisch et al.l 
(l20ri]l for a review of this work. The Sun is currently 
near the edge of the Local Interste llar Gloud (LIC). Ob¬ 
servations by iLinskv et al.l (j2000l l placed a lower limit 
on the neutral hydrogen number density of the LIC of 
uhi ^0.10 cm“^; this may be regarded as a firm lower 
limit for the total column density, as H i will be only one 
constituent of the LIC. The combination of more detailed 


models and observations find total hydro gen number den¬ 
sities ranging from 0.23 to 0.3 0 cm“^ (|Slavin fc Frischl 
120081 : iHeerikhuisen et al.l[201411 . We will thus adopt 0.2 
cm“^ as the LIC density for the remainder of this work. 

1.2. Effects of the ISM upon Planets 

The motion of planet-hosting stars through the ISM 
can have several effects upon the planets, due to the 
changing astrosphere size (which affects the cosmic ray 
flux experienced by the planet) as well as direct depo¬ 
sition of interstellar material onto the planetary atmo¬ 
sphere. We will briefly summarize these processes. 

The cosmic ray flux can have a number of ef¬ 
fects on the climate and habitability of planets 
(see IScherer et al.l 120061 for a review of these pro¬ 
cesses). Cosmic rays h ave been tentatively linked to 
cloud nucleation (e.g.. iSvensmark fc Friis- Chris tenseiil 
119971 : iMarsh fc Svensmarkll2000l : iKirkbv ^al.ll20lill. al- 

though this is controversial le.g.. iKernthaler et al.lll999l : 

iLaken fc Calogovi^l2013ll . lErlvkin fc Wolfendalel (120111 1 
surveyed the literature and concluded that the contribu¬ 
tion of cosmic rays upon the Earth’s tropospheric cloud 
cover is of order 1%. Lightning may also be affected 
by cosmic rays, as cosmic ray-induced atmosph eric ion¬ 
izatio n could affect the global electric circuit (jChronisI 
1200911 . Both clouds and lightning can have broader ef¬ 
fects on climate; clouds change the planet’s albedo, thus 
affecting surface temperatures, while lightning produces 
NOx compounds in the troposphere, which ar e precursors 
to th e formation of ozone, a greenhouse gas (|Allen et al.l 
|20I3). 

More extreme effects will occur if two of the Earth’s 
cosmic ray shields are lost simultaneously-that is, if 
the Earth’s magnetic field undergoes a polarity rever¬ 
sal when the solar system is traversing a d ense interstel¬ 
lar clou d. This scenario was considered bv I Pavlov et al.l 
(|2005all . During such a reversal, the overall strength of 
the Earth’s magnetic field decreases. The resulting in¬ 
crease in cosmic rays causes ionization of atmospheric 
molecules, which directly forms NOx compounds. Un¬ 
like in the troposp here, stratospheric oz one tends to be 
destroyed by NOx. iPavlov et al.l (|2005al) found that this 
could decrease ozone concentrations by 40% or more 
globally. For life-bearing planets, the cosmic ray flux 
has a direct bearing on mutation rates and further radi¬ 
ation damage, both directly through impacts of cosmic 
rays and secondary particles from cosmic ray air show¬ 
ers on living tissue, an d indirectly through reduction of 
the global ozone layer (IPavlov et al.ll200^ . resulting in 
higher doses of ultraviolet radiation. They found that, 
given the duration of such cloud crossings and the fre¬ 
quency of magnetic field reversals, a reversal should occur 
at least once during each cloud crossing. 

Interstellar dust can also affect the climate of planets. 
If sufficient amounts of dust are swept up into a planetary 
atmosphere, it can cause a “reverse greenhouse” effect, 
radiating efficientlY i n the infrared and cooling the planet 
(| Pavlov et aLll^ObbU . 

The most extreme possible effects will occur if the ISM 
density is sufficiently high that the astrosphere retreats 
inside the habitable zone . This possibility was inves¬ 
tigated by iSmith fc Scalol (|2009(1 , who referred to this 
occurrence as “descreening.” They found that for a Sun- 
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like star an ambient ISM density of ~600 cm“^ is nec¬ 
essary for descreening to occur. Clouds with such high 
densities are small, and they calculated a rate of these 
events of 1 — 10 Gyr“^ for sun-like stars, and less fre¬ 
quently for later-type stars due to th e greater proximity 
of the habitable zone to these stars. lYeghikvan fc Fahn 
(|2004fl found that if the Earth is directly exposed to the 
ISM, the increased flux of hydrogen could alter the at¬ 
mospheric chemistry, resulting in severe ozone depletion. 
This could also result in an ice age. 

1.3. Relating Observables to Astrospheres 

The Kepler mission searched for planets down to the 
size of the Earth using the transit method, by mak¬ 
ing nearly continuous photometric observations of nearly 
150,000 stars for more than four years and searching for 
the periodic flux decrements ass ociated with transiting 
planets (e.g.. iSeader et all 1201511 . The region that Ke¬ 
pler surveyed for planets (hereafter denoted the ‘‘‘‘Kepler 
field of view” or the “Kepler search volume”) is centered 
at an RA of 19'‘22™40'* and a Declination of 44°30'00", 
and fits within a square 13.9° on a side, although the 
filling factor of the Kepler CCDs w ithin this square is 
less than unity (jBorucki et al.l 1200811: the Kepler CCD s 
cover an area of 115 square degrees (ISeader et al.ll201^ . 
Kepler is the first mission capable of detecting poten¬ 
tially habitable planets-i.e., those with masses close to 
that of the Earth and orbiting far enough from their host 
stars that they could possess surface liquid water-where 
the effects of interactions with the ISM described above 
could have interesting consequences. Kepler has already 
found a number of sma ll planets in the habitable zone 
(e.g., iTorres et al.l 1201^ . We thus chose this region of 
the sky for a survey to determine the structure of the 
ISM, with the objective of constraining the interstellar 
environment for the known potentially habitable planets 
discovered by Kepler. 

The first to co nsider the ast rosphere sizes of exoplane¬ 
tary systems was iFrischl (1199311 . who calculated the astro- 
sphere size for a sample of 60 nearby G-type stars. At the 
time, however, none of these stars was known to possess 
planets. With the multitude of exoplanets now known, 
in particular the large number of candidate planets from 
Kepler^ including a number of confirmed habitable zone 
super-Earths, we can today make a much more specific 
study of the ambient ISM properties for known planetary 
systems. To our knowledge the present work is the first 
to address the ISM environments of known exoplanetary 
systems. In analogy to many recently-coined terms re¬ 
ferring to extrasolar equivalents to features of the solar 
system, we propose the term exoLISM to refer to the 
ISM in the vicinity of exoplanetary systems. In com¬ 
plementary work, Edelman et al. (2015, in prep) have 
searched for Ly a absorption from the hydrogen walls of 
astrospheres associated with nearby planet-hosting sys¬ 
tems in order to estim ate the mass loss rates (i.e., the 
stellar wind properties: iWood et al.ll2005ll of these stars. 
Most of the Kepler systems, however, are too distant 
for this kind of observations, as absorption from inter¬ 
vening interstellar gas overwhelms the weak astrospheric 
signal. Therefore, mapping the ISM in the Kepler search 
volume is currently the best option for constraining the 
astrosphere sizes for Kepler planets. 

We note that our techniques are only capable of prob¬ 


ing two of the four parameters most responsible for de¬ 
termining the size of an astrosphere (p/sm, vism, Psw, 
vsw, where “SW” denotes the stellar wind). We can 
measure the ISM density and the ISM velocity, the latter 
of which is necessary to determine the relative velocity 
between the ISM and the star. Our methods, however, 
are not capable of probing the transverse velocities, and 
hence the full 3-D space motions, of ISM clouds, and 
so we will have to estimate the full 3-D relative velocity 
from the difference in radial velocities between the star 
and the ISM. In this work we do not measure the other 
two parameters upon which the astrosphere size is most 
strongly dependent, namely the stellar wind velocity and 
density; these can be measured for nearby stars through 
astrospheric absorption, as described above. The astro- 
sphere size is also dependent upon the other properties 
of the ISM, such as the interstellar m agnetic field and 
ionization properties (|Zank et al.ll2ni3ll . but we are also 
unable to determine these parameters from our data. 

This work has a number of secondary applications in 
addition to estimation of astrosphere sizes. Eirst, the dis¬ 
tances to many of the Kepler candidate systems and even 
confirmed planetary systems are rather uncertain. De¬ 
termination of the locations and distances to ISM clouds 
within the Kepler field of view can help constrain the 
distances to these systems; if absorption from a given 
cloud is seen in the spectrum of a planet host star, then 
that star is at least as distant as that cloud. Conversely, 
if a star is located in the same part of the Kepler held 
of view as a cloud but does not show interstellar absorp¬ 
tion in its spectrum, then the star is at most as distant 
as that cloud. Second, this work will allow the prediction 
and modeling of interstellar lines in the spectra of Ke¬ 
pler targets, potentially an important source of spectral 
contamination for some applications. 

2. OBSERVATIONS, DATA REDUCTION, AND MODELING 

Targe ts were selected from the updated Hipparcos cat¬ 
alog of Ivan LeeuwenI (1200711 . We selected stars in the 
Kepler held of view with B — V < 0.41 in order to se¬ 
lect early-type stars, and concentrated on bright stars 
{V < 7) in order to minimize the necessary observation 
time; ^ 100 stars met these criteria. From this list we ob¬ 
served 17 targets; these consisted of the brightest stars 
in the sample, plus others selected to provide roughly 
even sampling over the full range of possible distances, 
from ^ 20 pc to ^ 1 kpc. The emphasis on very dis¬ 
tant stars resulted in a distribution that mostly sampled 
the low Galactic latitude side of the Kepler held at large 
distances. The properties of observed stars are listed in 
Table [TJ 

The data were obtained using the 2.7m Harlan J. 
Smith Telescope and Robert G. Tull TS21 Spectrograph 
(|Tull et al.lll995l l at McDonald Observatory on 2010 Au¬ 
gust 16—18 UT. The instrumental resolving power is 
R ^ 240,000. The weather varied from clear to light 
clouds. Two stars (HIP 96288 and 96693) were observed 
on both August 16 and 18 due to suboptimal signal-to- 
noise in the hrst observations. All other stars were ob¬ 
served on only one night. Our primary spectroscopic 
target was the Na i D1/D2 doublet, at AA 5896, 5890 
A. We also observed the K i line at 7699 A. Due to the 
instrumental setup, however, we could not observe the 
other line of the K i doublet. 
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Our reduction and analysis pip eline is the sam e as 
used on previous studies, such as iRedfieldl (|2007ll and 
IWvman fc Redfieldl (|2013il . We performed the data re¬ 
duction using standard IRAF tasks. We converted the 
velocities to the local standard of rest (LSR) frame, 
using the values of the solar motion measured by 
iSchdnrich et al.l (|201Clll . All velocities quoted in this work 
are in the LSR frame, unless noted otherwise. After 
aperture extraction, we normalized the data by fitting 
a polynomial function using sigma-clipping to isolate the 
continuum. This process was straightforward due to the 
flat continua of the hot stars that we targeted. Addi¬ 
tionally, since we have observed both lines of the Na i 
doublet we can minimize the effects of systematic errors 
in the normalization by simultaneously fitting both lines. 

Next we removed the telluric lin es in the data. We 
follow ed the same methodology as IWvman fc Redfieldl 
(1201311. using t he fo rward modeling techniques of 
iLallement et al.l (|1993ll with a model telluric spectrum 
produced by th e Atmospheric Transmission program 
(|Grossm^ 1198911 . The region around both Na i lines 
suffers from pervasive telluric contamination; the deep¬ 
est telluric lines absorb at most ^ 40% of the signal in one 
pixel. For the D1 line, the region between —15 km s“^ 
and -|-25 km s“^ is largely free of telluric contamina¬ 
tion. For the D2 line, however, there are telluric fea¬ 
tures that are usually superposed upon the interstellar 
features. Again, fitting both the D1 and D2 lines simul¬ 
taneously helps to minimize the effects of any imperfec¬ 
tions in the telluric subtraction. There is no significant 
telluric contamination within ^ 80 km s“^ of the K i 
line center. For some targets we needed to perform a 
second normalization after the telluric subtraction, us¬ 
ing a lower-order polynomial. 

We produced model spectral lines using a Voigt line 
profile; for the Na i lines we took into account the hy- 
perfine structure of each line. We then convolved this 
model of the actual line profile with the spectrograph 
line spread function to produce a model of the observed 
spectral lines. We performed fits to the D1 and D2 lines 
individually and simultaneously. We determined the fi¬ 
nal line parameters by combining the values from these 
three fits. For each line component, the free parame¬ 
ters are the velocity u, the Doppler parameter 6, and 
the column density N. We obtained error estimates us¬ 
ing a Monte Carlo routine after the fits were performed. 
We used the F-test to determine the best-fitting number 
of line components. For cases where no interstellar ab¬ 
sorption is detected, we calculated the upper limits on 
the column density from the root-mean-squared (RMS) 
scatter of the continuum in the region around the lines. 
For Na i, we use the D2 line to calculate the limits due 
to its higher oscillator strength. 

Some modifications of this process were necessary for 
HIP 96441, our coolest star (spectral type F4V), due 
to the presence of strong stellar Na i and K i lines in 
the spectrum. We produced an empirical model of the 
stellar lines by smoothing the spectrum and reflecting 
the resulting smoothed spectrum across the line center. 
These model lines were then subtracted from the data. 
This process removes the broad stellar lines while leaving 
any narrow interstellar lines intact. No interstellar ab¬ 
sorption, however, is detected for this target, the nearest 
in our sample {d = 18.34 pc). This is expected, as this 


target is well inside th e boundaries of the Local Bubble 
(|Lallement et al.ll2003ll . 

3. ANALYSIS 

We detect Na i towards 13 of the 17 targets (totaling 
39 line components), and K i towards 5 of the 17 (8 
total line components). All K i components correspond 
to a detected Na i component to within 1.2 km s“^. We 
show the data and fits in Fig. [T] for Na i and Fig. 
for K I. A number of the components, especially for 
the more distant targets, are saturated and/or blended 
with other components. The number of components and 
the total column density generally increase with distance. 
Detected Na i components cover a wide range of LSR 
velocities, from —14 km s“^ to -1-22 km s“^; all except 
three components have —10 km s“^ < v < -1-13 km s“^. 

The fit properties for all components are given in Ta¬ 
ble [2] for Na i and [3] for K i. The distributions of the 
velocities, Doppler parameters, and column densities for 
the Na i components are shown in Fig. [3^, b, and c, re¬ 
spectively. The velocities fall into a number of clumps, 
which we use to separate t he de tected Na i components 
into different clouds (see The majority of the 

detected ISM components have relatively narrow lines 
(6 < 2 km s“^), although three components have b ^ 4.5 
km s“^; the distribution shows two peaks, at ^ 0.5 km 
s“^ and ~ 2 km s“^. The distribution of Na i col¬ 
umn densities peaks at log(VNa i/cm“^) ^ 11.5, with 
a high column density tail of the distribution reaching 
to log(VNa i/cm“^) ^ 14. There does not appear to 
be an obvious correlation between component velocity 
and either Doppler parameter or column density. For 
comparison, the distributions of b and column density 
that we find for Na i ar e similar to those found by 
IWvman fc Redfieldl (|2013f) in the direction of the past 
solar trajectory; however, they measured a handful of 
components with higher Doppler parameters than we do 
(up to 5 ^ 8 km s“^), while we find some components 
with higher column densities than they did (their highest 
are log(VNa i/cm”^) ^ 12.3). 

Fig. [IK and b shows the total column density of Na i 
and K i, respectively, as a function of distance. As is ex¬ 
pected, the column density of both tracers generally in¬ 
creases with distance. There is, however, a large amount 
of scatter for Na i: ^ 1 dex for d < 200 pc, and greater 
than 2 dex at distances of ^ 700 pc. This is indicative 
of structure in the ISM. As there are large differences 
in the column densities towards targets with similar dis¬ 
tances separated by only a few degrees, there must be 
different distributions of Na i (i.e., different structures) 
along these lines of sight. We have several pairs of lines 
of sight separated by ~ 1° (corresponding to a physical 
separation of 1.7 pc at a distance of 100 pc or 12 pc at 
a distance of 700 pc), and many lines of sight showing 
different column densities separated by up to 14° (corre¬ 
sponding to 25 pc at a distance of 100 pc, or 175 pc at a 
distance of 700 pc). Therefore, we are sensitive to ISM 
structure transverse to the line of sight on scales of a few 
to tens of parsecs. __ 

Apparent optical depth (jSavage fc Sembachlll991h dia¬ 
grams for Na I and K i are shown in Fig. [31 see the figure 
caption for a full description of this plot. These plots 
do indicate some continuity of cloud structure; however, 
due to the relatively large size of the Kepler field and the 
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Table 1 
Target Sample 


HIP 

(1) 

HD 

(2) 

Name 

(3) 

RA (J2000.0) 
(4) 

Dec (J2000.0) 
(5) 

Spec. Type 
(6) 

V 

(7) 

d (pc) 

(8) 

Exp. time (s) 

(9) 

SNR 

(10) 

96441 

185395 

S Cyg 

19 

36 26.53 

+50 13 15.96 

F4V 

4.5 

18.34 + 0.05 

900 

130 

93408 

177196 

16 Lyr 

19 

01 26.37 

+46 56 05.32 

A7V 

5.0 

37.4 + 0.2 

1200 

110 

97165 

186882 

5 Cyg 

19 

44 58.48 

+45 07 50.92 

B9.5IV+ 

2.9 

50.6 + 1.2 

360 

160 

92822 

175824 


18 

54 47.12 

+48 51 33.84 

F3HI 

5.8 

55.3 + 0.8 

1200 

66 

97700 

188074 


19 

51 19.38 

+47 22 38.05 

F2V 

6.2 

62.5 + 1.6 

1200 

53 

96286 

184977 


19 

34 39.86 

+48 09 52.38 

A9V 

6.8 

86.5 + 2.8 

1200 

47 

96195 

184603 


19 

33 36.44 

+38 45 43.10 

A3Vn 

6.6 

129 + 8 

1200^^ 

51 

96288 

184875 


19 

34 41.26 

+42 24 45.04 

A2V 

5.3 

177 + 6 

1500‘> 

54 

96693 

185872 

14 Cyg 

19 

39 26.49 

+42 49 05.81 

B9IH 

5.4 

200 + 8 

1500‘> 

72 

93210 

176582 

V545 Lyr 

18 

59 12.29 

+39 13 02.36 

B5IV 

6.4 

292 + 26 

1500 

82 

98194 

189178 


19 

57 13.87 

+40 22 04.17 

B5V 

5.5 

340 + 22 

1200 

82 

96491 

185330 


19 

36 56.65 

+38 23 01.77 

B5II-IH 

6.5 

361 + 40 

1200 

53 

94481 

180163 

V Lyr 

19 

13 45.49 

+39 08 45.48 

B2.5IV 

4.4 

426 + 24 

900 

130 

97845 

188439 

V819 Cyg 

19 

53 01.25 

+47 48 27.79 

BO.SIIIn 

6.3 

503 + 71 

1200 

38 

97634 

187879 

V380 Cyg 

19 

50 37.33 

+40 35 59.14 

B1HI+ 

5.7 

649 + 101 

1200 

68 

95673 

183362 

V558 Lyr 

19 

27 36.40 

+37 56 28.31 

B3Ve 

6.3 

725 + 152 

1500 

43 

97757 

188209 


19 

51 59.07 

+47 01 38.42 

09.5Ia 

5.6 

1100 + 266 

1200 

72 


Note. — (1) Star ID in the Hipparcos catalog, in order of increasing distance. (2) Star ID in the Henry Draper catalog. 
(3) Bayer, Flamsteed, or variable star designation. (4) Right ascension, equinox J2000.0, in hours, minutes, seconds. (5) 
Declination, equinox J2000.0, in degrees, minutes, seconds. ( 6) Spectral type. (7) V^-band magnitude. (8) Distance in parsecs, 
as determined from Hipparcos parallaxes llvan Leeuwenlf2007n . (9) Total exposure time for the star, divided into 2-4 exposures 
per star. (10) Signal-to-noise ratio per pixel in the region of the Na I lines. Data for columns 1-8 from the SIMBAD database 
(http: / / simbad.u-strasbg.fr/simbad/). 

^ A cosmic ray hit on the Na I D1 line ruined one exposure; effective exposure length for analysis of the D1 line was 600 seconds. 
^ These targets were observed on two nights due to low signal-to-noise in the first set of observations due to thin clouds. All 
other targets were observed on only one night. 


aforementioned small-scale structure of the ISM, not all 
components which we infer to belong to the same cloud 
(see 1 13.11) are actually connected in this plot. 

The distribution of the detected Na i and K i compo¬ 
nents in velocity space and on the plane of the sky in 
relation to the Kepler field is shown in Fig. [6] As men¬ 
tioned above, most of our targets, especially those at 
large distances, lie on the side of the Kepler field closest 
to the Galactic plane (towards the lower left hand side 
of Fig. [5]). Thus, we only have limited information about 
the ISM over the higher Galactic latitude side of the Ke¬ 
pler field of view. Additionally, all five of the stars for 
which we detect K i absorption are located in the corner 
of the Kepler field closest to the Galactic plane. 

3.1. Identification of ISM Clouds 

Due to our relatively sparse spatial coverage we con¬ 
centrate on identifying ISM clouds based on the line ve¬ 
locity centroids. We assume that clouds move as approx¬ 
imately solid bodies, such that if we see two stars with 
absorption at the same LSR velocity at opposite sides of 
the field of view, they would result from the same cloud. 
In principle there could always be more distant absorp¬ 
tion components which happen to share the same veloc¬ 
ity as an actual cloud, and would thus be misidentified 
as belonging to the nearer cloud. Given our sparse spa¬ 
tial coverage we cannot positively identify such contam¬ 
ination, so we will proceed with the above assumptions 
with the caveat that such velocity interlopers may be a 
problem on longer sightlines. An additional complica¬ 
tion is caused by Galactic rotation. At a distance of 426 
pc (beyond which we only have four targets), objects at 
rest with respect to Galactic rotation towards our high¬ 
est and lowest Galactic latitude targets will show radial 
velocities d iffering by ~ 2.5 km s“^ (using Eqn. 5.170 of 
iLandflOO^ . This velocity differential is large enough to 


potentially cause confusion regarding the cloud identifi¬ 
cations. Most components, however, would show veloc¬ 
ity shifts due to Galactic rotation of < 2.5 km s“^ due 
to smaller distances and less widely separated sightlines, 
and so we again proceed with the caveat that radial ve¬ 
locity differentials due to Galactic rotation could be an 
additional source of confusion for our cloud identifica¬ 
tions. In order to set distance limits on the clouds, we 
define the lack of absorption at a given velocity as the 
lack of presence of a cloud. There could still be absorp¬ 
tion that is too weak for us to detect; typically, we are 
sensitive to column densities of log(A^/cm“^) > 10.7 for 
Na I and log(A^/cm“^) > 11.0 for K i. Lack of detec¬ 
tion of absorption could be due to either low densities, 
resulting in a column density of Na or K below our de¬ 
tection limits, or the presence of Na or K but in higher 
ionization states (Na ii, K ii, etc.), indicating high tem¬ 
peratures or the presence of ionizing radiation. Thus, we 
are only sensitive to a certain type of cloud: relatively 
cold (T < 1000 K), predominantly neutral clouds (e.g., 
iHobbsl 119781: iLallement et al.l 1200311 . Gonveniently, this 
is the type of cloud that is capable of significantly com¬ 
pressing astrospheres. Previous works to map the ISM 
using absorption l ine measurements hav e made similar 
assumptions fe.g.. ILallement et al.ll2003ll . 

As is apparent in Fig. [3^, there are six distinct clumps 
in the distribution: between —10 and —6, —5 and —1, 0 
and +2, -|-3 and +6, +6 and +9, and +9 and +13 km 
s“^. Each of these we identify as clouds due to their co¬ 
incidence and isolation in the velocity space, and their 
spacial continuity (see Eigs.|n]and[7]). We therefore iden¬ 
tify eleven ISM clouds in the Kepler field of view, five 
of which are only seen in absorption in the spectrum of 
one star. We attribute the other thirty-four detected 
Na I components to six other clouds, within the velocity 
ranges defined above. 
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No I D2 No I D1 




- 20-10 0 10 20 30 - 20-10 0 10 20 30 
Vlsr (km s“') Vlsr (km s“') 


Figure 1. Normalized spectra (thin lines) and fits (thick lines) at the Na i D2 (left) and D1 (right) lines for all stars. Spectra are shown 
in order of increasing stellar distance from top to bottom, and each spectrum is offset by unity for clarity. Vertical red hatch marks denote 
the central velocity of each component. Each spectrum is labeled with the stellar HIP catalog number. Spikes and dips in the spectra which 
are not included in the fits are largely due to imperfectly-removed telluric lines. These can easily be distinguished from Na i absorption 
because there is no corresponding feature in the other Na i line. 


We will now discuss the various identified clouds and 
single components, in order from most negative to most 
positive velocities. The bluest component is seen in the 
spectrum of HIP 97634 and is detected in both Na i and 
K I, with LSR velocities of —13.53 km s“^ (Na i) and 
— 13.43 km s“^ (K i). As there are no components de¬ 
tected near this velocity for other stars, we can only set 
an upper limit on the distance to this cloud of 649 pc, the 
distance to HIP 97634. The spatial extent is largely un¬ 
constrained, but it is not seen in the spectra of the only 
two more distant stars, HIP 95673 and 97757, limiting 
its extent towards the northeast and southeast sides of 
the Kepler field of view (the upper and lower left portion 


of Fig. m respectively). 

Next are the three Na i components between —10 and 
—6 km s“^ (seen in the spectra of HIP 97634, 97757, and 
98194), all located near the eastern corner of the field 
of view. Due to the spatial coincidence of these compo¬ 
nents, we identify this as “Cloud I.” The upper left panel 
of Fig. [7]shows the spatial distribution of these stars. We 
do not have any targets located in between these three 
stars, but no absorption in this velocity range is seen 
in the spectra of HIP 97700, located only 22’ from HIP 
97757. We thus consider the distance to HIP 97700, 62.5 
pc, to be the lower limit on the distance to this cloud, 
while the distance to the closest star with detected ab- 
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Table 2 

Na I Fit Properties 


HIP 

Comp. 

t'LSR (km s 1) 

b (km s ^) 

log(7VNa I/cm 2) 

Cloud 

( 1 ) 

( 2 ) 

(3) 

(4) 

(5) 

( 6 ) 

96441 

0 



< 10.7 


93408 

0 



< 10.5 


97165 

1 

-3.120 ±0.010 

0.268 ±0.047 


II 

92822 

1 

-3.96 ±0.32 

1.32 ±0.24 

ii.4tUO_Q Q24 

II 

97700 

0 



< 10.7 


96286 

0 



< 10.9 


96195 

1 

1.42 ±0.35 

0.231 ±0.052 

1U.OOO_Q QYQ 

III 

96288 

1 

0.89 ±0.11 

0.578 ± 0.029 

11.4296 ±0.0061 

III 

96693 

1 

-4.04 ±0.53 

2.9 ± 1.1 

10.789 ±0.019 

II 


2 

1.24 ±0.13 

0.64 ± 0.46 

lU.y fO_Q Q;l^g 

III 

93210 

1 

-4.19 ±0.20 

1.60 ± 0.10 

Q53 

II 


2 

0.867 ± 0.092 

0.76 ±0.23 

'^^-0.018 

III 


3 

3.65 ±0.16 

0.55 ±0.13 

li.DD4_Q Q24 

IV 

98194 

1 

-7.00 ± 0.49 

4.46 ±0.32 

10.909lg?gf 

I 


2 

-2.538 ±0.058 

1.35 ±0.10 

11.890t°;°20 

II 


3 

1.218 ±0.024 

0.99 ±0.13 

12.496l°;°g 

III 


4 

5.369 ±0.051 

1.11 ± 0.11 

-LO.40_q 43 

IV 

96491 

1 

-4.16 ±0.50 

1.65 ±0.22 

11 621 

li.DZl_Q Q45 

II 


2 

0.81 ±0.19 

1.47 ±0.19 

11.504 ±0.015 

III 

94481 

1 

-4.94 ±0.27 

2.45 ±0.77 

ii.iU_0 92 

II 


2 

1.079 ±0.022 

1.92 ±0.37 

11 097 + 0 .067 
LL.o^i _o.080 

III 


3 

8.194 ±0.087 

0.94 ±0.12 

11.902 ±0.014 

V 


4 

11.281 ±0.038 

0.298 ± 0.067 

J.X.‘±Z,y_Q Qgg 

VI 

97845 

1 

-3.74 ±0.34 

2.38 ±0.23 

iZ.UDU_Q 035 

II 


2 

3.98 ±0.15 

1.86 ±0.54 

12 39+0-3® 

xz.oy_Q 

IV 


3 

11.588 ±0.095 

1.90 ±0.17 

11.7897 ±0.0017 

VI 


4 

15.75 ±0.14 

0.44 ± 0.23 

11.76t°i® 


97634 

1 

-13.53 ±0.11 

1.882 ±0.072 

13.140t°;°?l 



2 

-6.639 ±0.065 

0.249 ± 0.099 

11.56t°iS 

I 


3 

-2.59 ±0.15 

0.50 ±0.33 

13.84t“;®^ 

II 


4 

1.09 ±0.32 

2.94 ±0.33 

12.218 ±0.056 

III 


5 

4.681 ±0.012 

0.378 ± 0.072 

12 039+0 066 
iz.uoy_Q 073 

IV 


6 

7.348 ± 0.018 

1.119 ±0.025 

11.997 ±0.016 

V 


7 

11.91 ±0.39 

1.47 ±0.41 

10.503 ± 0.037 

VI 

95673 

1 

-4.645 ± 0.056 

4.264 ± 0.052 

11.1861^32 



2 

-1.73 ±0.21 

0.44 ±0.11 


II 


3 

1.32 ±0.38 

0.56 ±0.37 

11 126^^'^^^ 
li.iZD_0 04g 

III 

97757 

1 

-8.63 ±0.13 

1.92 ±0.22 

12.53 ± 0.44 

I 


2 

-4.39 ±0.49 

0.83 ±0.25 

12.55 ± 0.79 

II 


3 

-2.13 ±0.95 

4.4 ±3.8 

11 075+0091 
ii.D/'0_o 115 



4 

4.939 ± 0.066 

0.866 ± 0.096 

12 07+0-26 
iZ.D/ _0 75 

IV 


5 

9.63 ±0.12 

1.71 ±0.13 

1 0 99Q + 0.0l3 
j.z,.zz,o_0_0i4 

VI 


6 

22.374 ± 0.095 

1.87 ±0.73 

10-786lg;°88 



Note. — (1) HIP catalog designation, in order of increasing distance. (2) 
Number of detected Na i component for that star (numbering from most nega¬ 
tive to most positive velocity), zero if no Na i is detected towards that star. (3) 
LSR velocity of component. (4) Doppler parameter of component. (5) Loga¬ 
rithmic Na I column density of component. (6) Cloud assignment of component 
(see text for details). 


sorption, HIP 98194 (340 pc), is the upper limit. Given 
the maximum separation between stars showing absorp¬ 
tion from Cloud I (6.7°), we can estimate the size of 
the cloud in the plane of the sky as 7.3 pc and 40 pc 
if located at the minimum and maximum distances, re¬ 
spectively. We note that this is a lower limit on the size, 
as the cloud may extend beyond the region probed by 
our survey. One Cloud I component, that towards HIP 
97757, is also detected in K i. 

Next is the largest clump in velocity space, eleven stars 
with absorption between —5 and —1 km s“^ (HIP 92822, 


93210, 94481, 95673, 96491, 96693, 97165, 97634, 97757, 
97845, and 98194), distributed across the field of view. 
There exist many stars with no absorption detected, with 
locations on the plane of the sky near or in between less 
distant stars with detected absorption. For instance, Na i 
absorption is not detected in the spectrum of HIP 97700 
or 96288, at distances of 62.5 and 177 pc, respectively, 
but is detected in the spectrum of HIP 97165, located al¬ 
most directly between HIP 97700 and 96288 in the plane 
of the sky, at a distance of 50.6 pc. 

We suggest two possibilities to explain the above obser- 
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Table 3 

K I Fit Properties 


HIP 

(1) 

Comp. 

(2) 

■ULSR (km s 1) 

(3) 

b (km s 

(4) 

log(ArK i/cm 

(5) 

An (km s 
(6) 

Cloud 

(7) 

96441 

0 



< 10.7 



93408 

0 



< 10.6 



97165 

1 

-3.214 ±0.051 

0.1008 ±0.0035 

10.206 ± 0.038 

0.014 ±0.052 

H 

92822 

0 



< 10.8 



97700 

0 



< 11.0 



96286 

0 



< 10.9 



96195 

0 



< 11.0 



96288 

0 



< 11.1 



96693 

0 



< 11.1 



93210 

0 



< 11.0 



98194 

1 

6.54 ± 0.11 

1.66 ±0.19 

10.847 ±0.032 

-1.17±0.12 

IV 

96491 

0 



< 11.0 



94481 

0 



< 10.6 



97845 

1 

4.96 ±0.069 

0.66 ±0.15 

11.377 ±0.066 

-0.97 ±0.17 

IV 

97634 

1 

-13.43 ±0.16 

1.54 ±0.29 

10.788 ±0.042 

0.10 ±0.19 



2 

-2.31 ±0.31 

0.94 ±0.60 

10.394 ± 0.078 

-0.28 ±0.34 

H 

95673 

0 



< 11.0 



97757 

1 

-8.743 ± 0.060 

0.202 ± 0.038 

12.39 ±0.32 

0.12 ± 0.14 

I 


2 

-4.90 ±0.23 

0.88 ±0.51 

10.420 ± 0.087 

0.51 ± 0.54 

H 


3 

4.633 ± 0.024 

0.393 ± 0.075 

11.59 ±0.11 

0.306 ± 0.070 

IV 


Note. — (1) HIP catalog designation, in order of increasing distance. (2) Number of detected 
K I component for that star (numbering from most negative to most positive velocity), zero 
if no K I is detected towards that star. Note that this is not the same as the component 
number of the corresponding Na I component listed in Table [2] (3) LSR velocity of component. 
(4) Doppler parameter of component. (5) Logarithmic K I column density of component. (6) 
Velocity difference between associated Na I and K I components, i.e., Av = UNa I — I- (7) 
Cloud assignment of component (see text for details). 


vations. First, this absorption may be due to one cloud, 
which contains a large amount of substructure. Alterna¬ 
tively, this absorption may be due to two or more smaller 
clouds which happen to share similar velocities. We do 
not have sufficient spatial or depth resolution to distin¬ 
guish between these possibilities, so for the time being 
we designate all the absorption between —5 and —1 km 
s“^ as “Cloud II,” with the caveat that this may not be 
one monolithic cloud. Three Cloud II components (those 
towards HIP 97165, HIP 97634, and HIP 97757) are also 
detected in K i. 

Due to the aforementioned substructure we cannot set 
a robust lower limit on the distance to this cloud, but the 
closest star with detected absorption is HIP 97165, at a 
distance of 50.6 pc. Thus, at least a part of this cloud 
is located at d < 50.6 pc, and is therefore the closest 
ISM structure found in this work and is located within 
the Local Bubble (the wall of the Local Bubble is lo- 
cated at a distance of ^ 100 — 150 pc in this direction; 
iLallement et al.l [20031) . The detection of this c loud at 
a distance of ^ 50 pc b y ILallement et aP (|2003l ) is dis¬ 
cussed in more detail in 93.31 Given the distance limits to 
the cloud and the separation of 13.9° between the most 
widely separated stars showing Cloud H absorption, we 
can estimate a size of the cloud of 12.5 pc if located at the 
maximum distance of 50.6 pc (smaller if located closer). 

We also detect very broad, shallow components in the 
spectra of HIP 97757 and HIP 95673 at velocities ap¬ 
propriate for Cloud H. In the case of HIP 97757 this 
component is not well separated from the other, stronger 
components seen in the spectrum, but its inclusion in the 
fit is supported by the F-test. These components have 
two of the largest Doppler parameters of any of our de¬ 
tected components (6 = 4.4 ± 3.8 km s“^ for HIP 97757 
and 4.264 ± 0.052 for HIP 95673). While these compo¬ 


nents could be produced by a single very hot or turbulent 
cloud, the very large Doppler parameters make this seem 
unlikely; these components may instead be caused by a 
number of narrower, unresolved absorption features, un¬ 
surprising as HIP 95673 and HIP 97757 are our two most 
distant targets, at 725 and 1100 pc, respectively. For 
these reasons we attribute the stronger, narrower Na i 
components in the spectra of these stars in the Cloud 
H velocity range to Cloud H, rather than these broad 
components. 

As is apparent from Fig. [3^, another, very narrow 
clump in the velocity distribution exists between 0 and 
-|-2 km s“^; the bluest and reddest components in this 
group are separated by only 0.6 km s“^ in velocity space. 
Na I absorption at these velocities is seen in the spectra 
of HIP 93210, 94481, 95673, 96195, 96288, 96491, 96693, 
97634, and 98194. The spatial distribution of these stars 
is presented in the upper right panel of Fig.jT] Again like 
Cloud I these stars form a spatially coherent group. We 
thus designate this as “Cloud HI.” All stars within the 
southeastern part of the Kepler field of view show absorp¬ 
tion at these velocities, and so we can only set an upper 
limit on the distance to this cloud of 129 pc, the distance 
to the nearest star with detected absorption, HIP 96195. 
The maximum separation between stars with Cloud HI 
detections is 11.2°, and so we can estimate a size limit of 
25.3 pc if the cloud is located at the maximum distance 
of 129 pc. None of these components are detected in K i. 

Absorption between -|-3 and -1-6 km s“^ is detected in 
the spectra of five stars (HIP 93210, 97634, 97757, 97845, 
and 98194). The spatial extent of these stars is shown in 
the lower left panel of Fig. [T] We refer to this as “Cloud 
IV.” Based upon the distances to the closest detection 
(HIP 93210) and the farthest non-detection in the region 
delineated by the detections (HIP 96693), we can set 
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Figure 2. Same as Fig.[T] except for the K I line. 

distance limits of 200 pc < d < 292 pc for Cloud IV. With 
a maximum separation of 13.0° between stars with Cloud 
IV detections, we can estimate sizes of 46.2 pc and 67.4 
pc if located at the minimum and maximum distances, 
respectively. Three of the five Cloud IV components, 
those towards HIP 97757, 97845, and 98194, are also 
detected in K i. 

Two stars (HIP 94481 and 97634) show Na i absorption 
between +7 and +9 km s“^. As can be seen in the lower 
middle panel of Fig. 0 these two stars both lie near the 
southern edge (bottom part of the figure) of the Kepler 
field of view, and there are no more distant stars with 
no absorption at these velocities located between these 
two stars. Thus, we tentatively label this as “Cloud V,” 
despite the detection towards only two stars; we note 
that in principle this could be two separate clouds which 
share the same velocity, but we will proceed under the 
assumption that this is a single cloud. The maximum 
distance to Cloud V is 426 pc, the distance to HIP 94481. 
We do not have any targets located directly between HIP 
94481 and 97634, and so we cannot directly set a lower 
limit on the distance to this cloud. However, the cloud 



-20 -10 0 10 20 30 

Vlsr (km s“') 



0 12 3 4 5 

b (km s“’) 



logf^No i/cm"^) 

Figure 3. a. Histogram showing the distribution of velocities of 
our detected Na I components. The velocity limits of our identified 
Clouds I through VI are also shown (see text for more details), b. 
Histogram showing the distribution of Doppler parameters of our 
detected Na I components, c. Histogram showing the distribution 
of column densities of our detected Na I components. Each com¬ 
ponent is colored according to its velocity in all three histograms. 

is unlikely to cover only the region of the sky directly in 
between the two stars. It cannot extend too far to the 
south, as no absorption is seen in the spectrum of the 
more distant HIP 95673 at these velocities. If we assume 
instead that HIP 94481 and 97634 are located near the 
southern edge of Cloud V, we can set a tentative lower 
limit of 200 pc through the non-detection of the cloud 
towards HIP 96288 and 96693. The separation of 7.2° 
between HIP 94481 and 97634 corresponds to cloud sizes 
of 25.3 pc at a distance of 200 pc and 53.8 pc at a distance 
of 426 pc. 

Four stars (HIP 94481, 97634, 97757, and 97845) dis¬ 
play Na I absorption between 4-9 and 4-12 km s“^. We 
attribute these four components to one cloud, “Cloud 
VI.” The spatial distribution of these stars is shown in 
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Figure 4. Observed Na i (a., top) and K i (b., bottom) column 
density towards each target star as a function of distance. Some 
error bars are smaller than the plot points. As expected, column 
density increases with distance, albeit with a large amount of scat¬ 
ter, indicative of small-scale structure in the ISM. 

the bottom right panel of Fig. [T] Based upon the closest 
detection and the nearest non-detection within this re¬ 
gion, we set distance limits of 200 pc < d < 426 pc. The 
maximum separation between stars with Cloud VI detec¬ 
tions is 8.77°, and so we can set limits on the transverse 
size of 30.7 pc and 65.3 pc if located at the minimum and 
maximum distances, respectively. No K i is detected for 
any of these components. 

Another single component is seen in the spectrum of 
HIP 97845 at a velocity of 15.75 km s“^. This cloud 
would lie at a maximum distance of 503 pc. Finally, the 
reddest component is detected in absorption towards HIP 
97757, at a velocity of 22.374 km s“^ and a maximum 
distance of 1100 pc. 

The ranges of parameters seen in each of our identified 
clouds with more than one observed sightline are listed 
in Table m 

3.2. Comparison of Na i and K i 

All eight detected K i components lie within 1.2 km 
s“^ of a detected Na i component. Five of the eight have 
velocities identical to the corresponding Na i component 
to within la. Almost without exception the Na i com¬ 
ponents with a corresponding K i component are those 
with the highest column densities for a given star. Three 
K I components are associated with Cloud IV, three with 
Cloud H, one with Cloud I, and one with a cloud detected 
along only one sightline. Additionally, the sightlines with 
detected K i absorption tend to be those located closer 
to the Galactic plane (off the lower left side of Fig. |6]). 

Fig. [8] compares the Na i and K i column densities for 


the targets where both are detected, along with the re- 
lation between the res pective column densities found by 
IWeltv fc Hobbsl (|2001h . The scatter around their rela¬ 
tion is large, though larger when measured by component 
(Fig. EK) than by sightline as a whole (Fig. |8b). We cal¬ 
culated 3cr upper limits on the K i column density that 
could be associated with each Na i component, under the 
assumption that &k i = ^Na i- As can be seen in Fig. El 
these upper limits are largely uninformative. 

Observations of two or more insterstellar species can be 
used to probe the temperature of and turbulence in the 
ISM through the Doppler parameter (line width). From 
iRedfield fc Linsk-(1 (|2004bll . the relationship is 

9lrT T 

6^ = —-b = 0.0166294-b (1) 

m A 

where k is Boltzmann’s constant, m is the mass of the 
given ion, ^ is the turbulent velocity, and A is the mass 
of the ion in atomic mass units; the last part of Eqn. [I] 
assumes velocities in km s“^. 

From Eqn. [1] it is apparent that, for a given i, there 
is a range of physically allowed values of 6k i, with the 
maximum corresponding to pure turbulant broadening 
(6k I = ^Na i) and the minimum corresponding to pure 
thermal broadening (6 k i = \/A-^a. i/Ak i6Na i)- These 
limits are shown in Fig. El 

Of our eight components with both detected Na i and 
K I, only one lies within the physically allowed region 
in Fig. Eland two others are consistent with this region 
to within Icr. The remaining five components are more 
discrepant. This could occur if the Na i and K i are 
in general not well-mixed, i.e., they do not share the 
same temperature and/or turbulent velocity distribution. 
Nonetheless, we note that all of our detected K i compo¬ 
nents lie within 1.2 km s“^ of a Na i component, indi¬ 
cating that some physical relationship exists between the 
two ions, e ven if they are not wel l-mixed. Alternatively, 
as noted bv lWeltv fc Hobbi (|2001ll . due to the differences 
in oscillator strength and abundance for Na i and K i, 
typically either the Na i lines will be saturated or the 
corresponding K i line will be very weak. Thus, our re¬ 
sults could also be explained if we have underestimated 
the systematic errors on the Doppler parameters due to 
the saturation of the Na i lines or the weakness of the 
K I lines. There is no correlation between the ratio of the 
6 values and the distance to the target (as might be ex¬ 
pected if blending of multiple clouds in longer sightlines 
results in corrupted line widths). There is also a large 
range in the ratio of the 6 values among the components 
attributed to each of Clouds H and IV. 

3.3. Comparison with Previous Work 

To our knowledge, only two of our target stars have 
been previously observed in ei t her N a i or K i. HIP 97165 
was observed by IWeltv et al.l (1199411 in Na i: our re sults 
are consistent with theirs. iChaffee fc Whitel (1198211 ob¬ 
served both lines of the K i doublet in absorption towards 
HIP 97757, albeit at lower resolution than the present 
work. They detected components at heliocentric veloc¬ 
ities of —25.0 and —11.4 km s“^, (or LSR velocities of 
— 10.1 and 3.5 km s“^), likely corresponding to our com¬ 
ponents at —8.473 and 4.663 km s“^, respectively. Their 
Doppler parameters (1.3/^J g km s“^ and 0.8 ± 0.1 km 
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Figure 5. Apparent optical depth plot for Na i (left) and K i (right). Each star is represented as a horizontal line; solid lines denote 
stars with detected Na i or K i, while dashed lines denote stars with no detected Na i or K i. The color scale denotes the column density 
per unit velocity in Na i or K i. At the location of each star this simply shows the observed absorption line profile, while elsewhere this 
distribution has been interpolated between the observed stars; this interpolation is necessary as we only have constraints at the locations 
of the horizontal lines. Note that these figures do not show the physical distribution of ISM clouds; the length of a contour in the vertical 
direction is not the physical size of the cloud along the line of sight. For example, a 1 pc thick cloud at a distance of 200 pc covering the 
entire Kepler field of view would manifest as a dark streak beginning at a distance of 200 pc and running unbroken to the top of the figure. 
Absorption shown at a given velocity and distance in this figure is caused by absorbing material that could be physically located at any 
smaller distance. The discontinuous colorscale at a given velocity is because of variations in the ISM across the Kepler field of view; a given 
star may or may not show absorption from any given cloud, dependin g upon its location in t he plane of th e sky. For Na i, the column 
density is also converted to total H column density using the relation of lWeltv &; HobFl l|200Ul : see also El 


s respectively) are somewhat larger than ours. Their 
measured logarithmic column densities (11.26±0.14 and 
ll.44lo.39j respectively), deviate from ours by 3.2cr and 
1.2(7, respectively. These discrepancies likely result from 
their lower spectral resolution. They did not detect our 
third component at a LSR velocity of —4.90 km s“^; 
with log Vk I = 10.420, it was likely below their detec- 
ti on limit. _ 

iLallement et al.l (j200,3| l presented maps of the neutral 
gas ISM density within ~250 pc derived from Na i equiv¬ 
alent width measurements and column densities. In the 
two maps in their work bracketing the Kepler field (see 
their Figs. 7 and 8), a small cloud is visible at a distance 
of ^ 50 pc. They identify this c loud w ith the Na i compo¬ 
nent identified bv lWeltv et al.l (| 19941) towards HIP 97165, 
and is thus the same as our Cloud II. They also identified 
absorption from this cloud towards HD 192640 (29 Cyg) 
and HD 193369 (36 Cyg), both located well outside the 
Kepler field of view, and quoted a heliocentric velocity 
of —19 km s“^ for these components (corresponding to 
LSR velocities of ~ —4 km s“^), in agreement with our 
values for Cloud H (between ^ — 5 and ^ — 1 km s“^). 
This suggests an angular size for the cloud of 19.3° and 
a ph ysical size of 17.7 pc at a distance of 50.6 pc. How¬ 
ever, ILallement et al.l (120031) also noted that HD 192640 


is a A Boo star, and so the Na i absorption could be 
ci rcumstellar in natu re. 

iWeltv et al.l (|1996D also observed HIP 97165, albeit in 
Ca II. They detected three absorption components, the 
weakest of which is coinci dent in velocity with our one 
detected Na i component. iRedfield fc Linsk-s^ (I2008D at¬ 
tributed this component to the Aql cloud, one of the 
LISM clouds close to the Sun. For other sightlines 
through the Aql cloud, they measured H i column den¬ 
sities of log(VH i) = 17.1 — 18.1. Assuming that the H 2 
fraction is negligible (reasonable for the modest density, 
hot conditions of the LISM), th e relation between Na i 
and H co lumn density found bv IWeltv fc Hobbl (|200H ) 
(see (13.41 for more detail) predicts that the Na i column 
(iensity towards HIP 97165 due to the Aql cloud should 
be log(VNa i) < 7.3, well below our dete ction limits. We 
note t hat the data used to derive the IWeltv fc Hobbsl 
(|2001l) relation do not extend below log(VH 1-1-H2) 19, 

and so this estimate should be considered very uncer¬ 
tain; nonetheless, it suggests that our detection of Na i 
towards HIP 97165 is not due to the Aql cloud, but rather 
due to th e more distant, higher (j olumn density cloud de¬ 
tected by ILallement et al.l (120031) , which is coincident in 
velocity with the Aql cloud. 
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22.37 


13.40 


4.42 


- 4.55 


- 13.53 



Figure 6. Map of the Kepler field, showing the observed stars and the Kepler CCD layout. Each star is scaled according to the inverse 
of its distance, with larger symbols denoting more distant stars. Each plot symbol is further divided into a number of slices, each denoting 
one ISM component detected in Na i and colored according to velocity using the scale at left, identical to that used in Fig. [S] white circles 
denote stars with no detected Na i. K i detections are denoted by small squares to the left of each star, colored according to the same 
velocity scale. The overlayed coordinate grid shows lines of constant RA and Dec, with a spacing of 20™ for the RA grid and 5° for the Dec 
grid. The Kepler field of view is centered at a = 19^22"^40^ S = +44°30^00^^ JBorucki et al.|[^08lh and each of the square CCD modules 
covers an area of 5 square degrees (http://kepler.nasa.gov/science/about/targetFieldOfView/). 


3.4. Estimating Volume Densities 

IWeltv fc Hobbsl (120011) found a relationship between 
Na I and total H column density, namely (from Fig. 18 
and associated text in that work), 

log Vh 1 +H 2 = 0.478 log iVNa I + 14.6 (2) 

which we adopt to calculate the total H column density 
from our Na i measur ements. We note, however, that 
IWeltv fc Hobbsl (120011) also found that even at a fixed 
Na I column density, the total H column density can vary 
by up to ~ 1 dex. Thus, our estimated column densities 
will not be particularly accurate. 

In order to produce physically motivated upper and 
lower limits on the number density of each cloud, we 
need to have physically motivated limits on the cloud 
size along the line of sight, which, when coupled with 
the observed column densities, will yield limits on the 
number densities. We emphasize that the limits on cloud 
size that we derive here are used only for computing the 
limits on number density, not for our later discussion 
of which planetary systems might he within each cloud. 
Observed ISM clouds typically are in the form of sheets 
or filaments (such as the interstellar cirrus discovered by 
IRAS; ILowet ahlH^ . or, sometimes, rounder, more 


“potato”-like shapes (e. g., the LIC that curren tly sur¬ 
rounds our solar system; [^dfieUi&IOsik^IinQS)- Thus, 

we can obtain a reasonable upper limit on the size of a 
given cloud (and a lower limit on the number density) 
by assuming that it is spherical, i.e., that the depth of 
the cloud in the radial direction is equal to its maximum 
extent in the tangential direction. We estimate this by 
calculating the separation between the two most distant 
stars for which absorption is detected for a given cloud, 
and computing the corresponding transverse size at the 
distance corresponding to the distance upper limit for 
that cloud. We also need a lower limit on the cloud size in 
order to ca l culate an upper limit on the number density. 
iPeek et ^ (|2011l) found that the cold, dense Local Leo 
Cold Cloud (LLCC), an island of cold neutral medium 
within the Local Bubble, has a width in the plane of 
the sky of ~ 0.25 — 0.54 pc; assuming that its thickness 
along the line of sight is similar to its width (i.e., that 
the cloud is tubular), they estimated a number density 
of 150-320 cm“^ for the cloud. Dense clouds and cores 
in s tar-formiiig reg ions can have thic knesses from 0.1 5 
pc (|Lee et al.ll201^ down to 0.01 pc (iWhite et al.ll201^ 
However, these clouds are very dense— iLee et al.l ( 20141) 
found densities of ^ 1—2 x 10® cm“®, while smaller clouds 
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Cloud IV, 200 pc < d < 292 pc Cloud V, 200 pc < d < 426 pc Cloud VI, 200 pc < d < 426 pc 



Figure 7. Same as Fig. [6] but showing each cloud identified in ED For each cloud only the stars more distant than the closest detection 
of that cloud are shown. Stars with absorption detected from that cloud are colored by velocity according to the same scale as in Fig. [6] 
while those with no absorption detected are shown in gray, in order to show the spatial continuity of each cloud. The coordinate grid is 
the same as in Fig. |6] 


will be even denser. Even with a thickness of 0.1 pc, our 
highest column density clouds would only have number 
densities of ^ 5 x 10^ cm“^, while if they have thick¬ 
nesses of 0.5 pc they would have densities of the same 
order of magnitude as the LLCC. We thus conclude that 
thicknesses of 0.5 pc are more plausible for the types 
of clouds that we have observed. Moreover, such dense 
clouds tend to be physically small, making it unlikely 
that we would see such a cloud stretching over a signifi¬ 
cant portion of the Kepler field of view. We thus adopt 
0.5 pc as a reasonable lower limit to the cloud thickness; 
the exact value for the lower limit is not overly impor¬ 
tant, as later we will only use the cloud density lower 
limits derived from the cloud size upper limits, not the 
density upper limits derived from the size lower limits, to 
estimate astrosphere sizes (see ^4.21) . Note that this says 
nothing about the orientation of the cloud with respect 
to the line of sight; such a cloud could be oriented at an 
angle to the line of sight, and so have considerable overall 
depth (i.e., difference between the distances to the clos¬ 
est and farthest elements of the cloud) even if any given 
line of sight only transverses the cloud for 0.5 pc. 

We then combine these size estimates with the max¬ 
imum and minimum measured column densities for a 
given cloud to compute reasonable upper and lower lim¬ 
its for the number density in each cloud. In most cases 
these limits are rather unconstrained, but are consistent 
with measured volume densities in the ISM (see Table 
[4] for the complete list of upper and lower limits). The 
upper limit for Clouds II is ^ 10^ cm“^, likely an over¬ 


estimate as this density would correspond t o a molecular 
cloud, and the survey of iDame et al.l (|200Ill did not find 
any significant CO emission in the direction of the Kepler 
field. 

For all six clouds, the minimum number densities that 
we derive are equal to or greater than that inferred for the 
Lie, the cloud currently surrounding our solar system. 
Thus, any planetary systems around solar-type stars that 
currently reside within these clouds can be expected to 
possess astrospheres smaller than the current heliosphere 
(for a similar velocity differential between the star and 
its exoLISM and similar interstellar magnetic field and 
ionization properties). 

4. DISCUSSION 

4.1. Placement of Planets in Clouds 

Many of the confirmed Kepler planets do not have re¬ 
liable published distances, due in no small part to the 
faintness of the host star population on average. A hand¬ 
ful, however, do have distances determined through a 
variety of methods. In Table [S] we summarize all con¬ 
firmed or validated planets in the Kepler field of view 
with known distances of less than 450 pc. We choose 
this distance limit as 426 pc is the upper limit for the 
distance to Clouds V and VI, the most distant of our 
clouds observed on more than one sightline. 

In addition to the Kepler planets, there is one known 
and one suggested planet from radial velocity surveys in 
the K epler field of view. 16 Cyg B b (iWittenmver et ^ 
1200711 is an msini = 1.68Mj planet on an ~ 800 day or- 
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Table 4 

Range of Na I Cloud Properties 


Cloud Name 
(1) 

too 

'•— 

«LSR (km s 1) 

(3) 

b (km s ^) 

(4) 

log(ArNa I/cm 2) 

(5) 

i/N-k I 
(6) 

nn (cm 

(7) 

Distance (pc) 
(8) 

I 

3 

-8.63-6.639 

0.249 -4.46 

10.909 - 12.53 

1.4 

0.51 - 260 

62.5 - 340 

II 

11 

-4.94-1.73 

0.268- 2.9 

10.789 - 13.84 

29 - 2800 

1.5 - 1100 

< 50.6 

III 

9 

0.81 - 1.42 

0.231 - 2.94 

10.583 - 12.496 


0.58 - 240 

< 129 

IV 

5 

3.65 - 5.369 

0.378 - 1.86 

11.664 - 13.45 

10 - 400 

0.73 - 690 

200 - 292 

V 

2 

7.348 - 8.194 

0.94- 1.119 

11.902 - 11.997 


1.2 - 140 

200 - 426 

VI 

4 

9.63 - 11.91 

0.298 - 1.90 

10.503 - 12.190 


0.21 - 170 

200 - 426 


Note. — (1) Designation of cloud (see EU- (2) Number of sightlines along which the cloud is detected in Na i. 
(3) Range of observed LSR velocities for this cloud. (4) Range of Doppler parameters for this cloud. (5) Range of Na i 
column densities for this cloud. (6) Range of Na i to K i column density ratios (i.e., abundance ratios) for this cloud. No 
detected K i co mpo nents were associated with Clouds I, III, or V. (7) Range of possible hydrogen number densities for 
this cloud (see ^3.411 . (8) Range of possible distances to the cloud. 


0. By component 




log(NK ,/cm“^) 

Figure 8. Comparison of the column density in Na i and K i 
for each component where Na i is detected (a., top) and for each 
sightline (b., bottom). Some error bars are smaller than the plot 
points. Upper limits in Na i or K i are shown as arrows. The line 
denotes the re lationship between the column densities of these two 
ions found by I Wei tv fc Hobbs! II2001I) . 

bit located at a distance of 21.21 pc (Ivan LeeuweDll200^ . 
This is sufficiently close that it likely res ides well within 
the Local Bubble (|Lallenient et al.l 1200,11 1 . and likely ex¬ 
periences a modest density interstellar environment sim¬ 
ilar to that of the Sun. 

The closest of our target stars, HIP 96441 (a.k.a. 
0 Cyg), has been suggested to harbor one or more plan¬ 
ets, but an unusual correlation between the bisector ve¬ 
locity span and the radial velocities m akes this inter¬ 
pretation ambiguous (|Desort et al.ll200^ . Located at a 
distance of 18.34 pc, 9 Cyg is likely also not surrounded 
by any of the ISM clouds detected in this sample, by 
the same reasoning used for 16 Cyg. Indeed, we do not 
detect any absorption towards this target. 



Figure 9. Relationship between the Doppler parameters for Na i 
and K i for our data (points) and the theoretically expected range 
(between the lines correspondings to purely turbulent broadening, 
top, and purely thermal broadening, bottom). Some error bars are 
smaller than the plot symbols. 

Most of the nearby Kepler planets are on orbits close to 
their stars and thus are not in the habitable zone. There 
are, however, a few interesting systems in this sample. 
Most interesting for our purposes are Kepler-186, host 
of the first confirmed appr oximately Earth-size habitable 
zone planet, Kepler-186 f (I Quint ana et al.ir2ni4[l . Kepler- 
22, host of the first confir med Kepler planet i n the hab¬ 
itable zone, Kepler-22 b (iBorucki et al.l[2Mll . and sev¬ 
eral systems with small habitable zo ne planets statisti¬ 
cally validated bv iTorres et al.l (1201511 : Kepler-438, 296, 
440, 441, 442, 437, and KOI-4427. We choose to in¬ 
clude KOI-4427 in our sample, even though the planet 
candidate KOI-4427.01 has only been statistically vali¬ 
dated through the exclusion of false positive scenarios 
to a confi dence level o f 99.2% , unlike the other sys¬ 
tems from iTorres et al.l (1201511 . which have been vali¬ 
dated to a confidence level of 99.5% or higher. Other sys¬ 
tems include Kepler-16, host to the first confir med tran¬ 
siting circumbinary planet (iDovle et al.ll201lll : Kepler- 
20, a system with five tra nsiting super-earths and nep- 
tunes l|Gautier et al.l[20T^ : Kepler-42, a system of three 
sub-Ea rth-size planets i n ver y close orbits around an M 
dwarf (|Muirhead et al.l [201211 : Kepler-37, host to three 
planets, including the smallest known transiting planet, 
with a radius only sl ightly larger than that of the Moon 
(|Barclav et al.l 1^01311 : Kepler-421, host to the longest- 
period confirmed transiting planet, w ith an orbital pe¬ 
riod of 704 days (|Kipping et al.l [201^ : and Kepler-444, 
an ~11 Gyr old thick disk star with five transiting plan- 
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ets. In total, we consider 31 systems. 

For each of these systems, we considered their distances 
and locations in the Kepler field of view relative to our 
identified clouds to determine whether any of these plan¬ 
etary systems could be located within one of these clouds. 
We find that 13 of these systems could lie within one of 
our identified clouds. Many of the Kepler systems, how¬ 
ever, are located in the central or northwestern parts of 
the Kepler held of view, where we do not have any tar¬ 
get stars at distances similar to these exoplanet hosts. 
Thus, our current sampling is insufficient to determine 
whether these systems could be located inside our ISM 
clouds. See Table [5] for a summary. 


4.2. Estimating Astrosphere Sizes 

Given this discussion regarding the placement of Ke¬ 
pler systems within our ISM clouds, we wish to estimate 
plausible astrosphere sizes. The size of an astrosphere, 
however, depends not only upon the exoLISM density 
and velocity but also on the parameters of the outhow- 
ing stellar wind. The stellar wind parameters are very 
difficult to observe for s o lar-ty pe stars, due to the small 
mass flux. IWood et al.l (1200511 . however, measured the 
mass-loss rate M for several nearby stars and provided 
a scaling relation between M and the x-ray hux Fx (for 
Fx < 8 X 10^ ergs cm“^ s“^). If we had x-ray huxes for 
the Kepler planet-host stars, we could thus infer M. 

We searched for x-ray emission from these stars by 
querying the NEXXUS 2 databas^l at the position of 
each of the host stars listed in Table [S] Unfortunately, 
none of these stars has been detected in x-rays. It would 
thus be helpful to obtain x-ray fluxes for Kepler planet- 
host stars in order to perform a more rigorous version of 
this analysis. Instead, we will use scaling relations based 
on the stellar mass and age in order to estimate the as¬ 
trosphere sizes. Ideally, the astrosphere sizes would be 
modele d using a detailed in ultifluid model, such as that 
used bv IMiiller et al.l (|2006f) to model the response of the 
heliosphere to differing LISM conditions; however, given 
the current uncertainty in the parameters of exoplanet 
host stars and their ambient ISM conditions, such a de¬ 
tailed a pproach is not warrant ed at this time. 

From iSmith fc Scalol (|2009ll . the size of a pressure- 
supported astrosphere moving supersonically and su- 
peralfvenically through the ISM, so that the dominant 
source of external pressure is ram pressure, is 


ro \n ) V 


( 3 ) 


where is the astrosphere radius, ro is an arbitrary 
reference radius within the astrosphere, uq and vq are the 
stellar wind number density and velocity, respectively, 
at ro, n is the IS M number densi t y, and V is the ISM 
streaming velocity. ISmith fc Scalol (120091 1 parameterized 
the dependence of the stellar wind on the age and mass 
of the star in a function 







4.5Gyr 


( 4 ) 


where M* and f* are the stellar mass and age, respec- 

http://www.hs.uni-hamburg.de/DE/For/Gal/Xgroup/nexxus/nexx 


tively. Thus, they had 


ro \ n J 


1/2 


V 


( 5 ) 


where ro, no, and Vg are now appropriate values for the 
Sun. We can write the solar wind density no in terms 
of the solar mass-loss rate, and thus obtain an absolute 
astrosphere size 




( 6 ) 


where m is the mean mass of a solar wind atom (which 
we assume to be equal to the mass of a hydrogen atom). 
This assumes that vq is constant over mass and age of 
the stars, which is a typical assumption o f previous works 
(ISmith fc Scald l^009[ IWood et al.l 1200511 . As discussed 
bv I WoodT i 200411 . this is probably a not unreasonable as¬ 
sumption given that the surface gravity of most solar- 
type main sequence stars is similar, but the wind velocity 
may be higher fo r rapid ly rotating stars. 

ISmith fc Scalol (1200911 noted that the yalues of the 
power law indi ces a and S are very uncertain. Us¬ 
ing d ata from IWood et al.l (|2005ll and IPenz fc Micelal 
(1200811 they calculated /? = —2.33 and /3 = —1.8, re- 
specti yely. We therefore ad opt /3 = — 2 for our calcula¬ 
tions. ISmith fc Scald (1200911 also noted that the yalue of 
a is e yen more uncertain; obseryations by IWood et al.l 
(j2005ll suggested that lower-mass and typically more ac- 
tiye stars haye stronger winds (i.e., a < 0). Once the 
x-ray flux increases beyond ^ 10® erg cm“^ s“^, how- 
eyer, a dramat ic drop in the mas s loss rate is obseryed 
for these stars IWood et al.1 (|2014l l. resulting in a much 
weaker wind (i.e., a > 0). With a lack of concrete in¬ 
formation, we adopt a = 0. For comparison we also 
calculated astrosphere sizes for a = —1,1; for sun-like 
stars the differences are negligible, but for the lowest- 
mass stars this results in an uncertainty of a factor of 
10 on the astrosphere size. This highlights the need for 
reliable x-ray flux measurements of planetary host stars, 
particularly for the lowest mass stars. 

For each star, we adopt parameters for the literature 
for M* and G. For the yelocity V of the star with respect 
to the ISM, we obtain the radial yelocity of the star from 
the literature, conyert this yelocity to the LSR frame, 
compute Av = |n* — 'CcioudI, and finally assume that V = 
VSAv, as the full three-dimensional space yelocity should 
be, on average, i/3 larger than the line-of-sight y elocity. 

In order to calibrate the ISmith fc Scald (1200911 model, 
we te sted it agains t the 2 7 multifluid heliosphere mod¬ 
els of lMfiller et al.l (1200611 by taking the releyant initial 
conditions for each of thes e models (titsm, utsm) and us¬ 
ing the m as inputs for the Smith fc Scffia (1200911 model. 
For the iMiiller et al.l (|2006 ) Model I we oyerpredict the 
heliosphere radius by a factor of 3.6; howeyer, this is 
unsurprising as the initial conditions for Model 1 are a 
yery hot, rarefied ISM, where the ISM p ressure is dom¬ 
inated by thermal pressure, whereas the ISmith fc Scalol 
(|2009ll model assumes that ram pressure is the domi¬ 
nant source of external pressure. For the other 26 mod¬ 
els, the disc repancies betwe e n the two models are much 
sm a.ller: the ISmith fc Scalol (|2009(1 model underpredicts 
^glfe Wl filler et alF i 200611 astrosphere radii by a factor of 
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1.14 to 2.18, with a mean of 1.54 and a standard de¬ 
viation of 0.24. We thus choose to m ake an empirical 
correction to the lSmith fc Scalol ( 2009 ) model to allow it 
to better match thelMiille r et al.l (l2nn(^ model results, by 
multiplying the iSmith fc Scald (1200911 astrosphere radi i 
by a factor of 1.54. This corrected ISmith fc Scalol (2009) 
model reproduces Models 2-27 of iMiiller et al.l (2006) 
with a standard deviation of 0.16. 

We also tested this empirically corrected model on 
the heliosphere , usin g the density of the LlC from 
iSlavin fc Frischl (j2008D an d the relative velocity b etween 
the Sun and the LlC from iMcComas et al.l (j2012t) . This 
yielded a heliosphere size of 154 AU, rather larger than 
the heliopause distance of 121 AU measured by Voyager 
1. Our calculations are for the heliopause distance in the 
upwind direction, and since Voyager 1 is not moving in 
this direction, the actual heliopause size with which we 
should compare our model will be < 121 AU. This sug¬ 
gests that our model is accurate to no better than 30%. 
Given the uncertainty in the hydrogen densities calcu¬ 
lated from the relation of iWeltv fc Hobb^ (120011) . how¬ 
ever, this level of inaccuracy in the astrosphere model 
will be only a small contribution to our error budget. 

We perform this calculation for each of the Kepler plan¬ 
ets which we find could be located within an ISM cloud. 
As noted earlier, we consider our upper limits on hydro¬ 
gen number density for most of our clouds to be much 
higher than reasonable values, and so for these calcu¬ 
lations we use our lower limits on uh for each cloud, 
resulting in upper limits on astrosphere size. These esti¬ 
mates are presented in Table[5] We were unable to locate 
measurements of the absolute radial velocities of HAT-P- 
11 (Kepler-3) and Kepler-45 in the literature, and so do 
not estimate astrosphere sizes for these systems. We also 
note that the value of /3 = — 2 gives a strong dependence 
on the stellar age; as accurately measuring stellar ages 
is notoriously difficult, these estimates should be treated 
with caution. For example, there are two s omewhat con- 
flicting ages in the literature for Kepler-20. iGautier et al.l 
() 2012 t) obtained an age of 8 . 8^2 7 Gyr from an isochrone 
analysis using the stel lar mean density m e asure d using 
the transit lightcurve. iWalkowicz fc Bas^ (I2013D calcu¬ 
lated an age of 4.0 Gyr using gyrochronology, where they 
measured the stellar rotation period using spot modula¬ 
tion in the Kepler lightcurve and used this and the rela¬ 
tionship between stellar age and rotation rate to mea- 
sure the stellar ag e . We adopt the isochrone age, as 
IWalkowicz fc Basril (|2013[ ) noted that the gyrochronolog- 
ical age-rotational period relations are significantly un¬ 
certain for stars of similar or greater age than the Sun; 
these relations are typically calibrated using open clus¬ 
ters of known age, and old open clusters are very rare. 
We note that the gyrochronological age for Kepler-20 is 
still compatible with the isochrone age to within 2a. Us¬ 
ing the gyrochronological age rather than the isochrone 
age gives astrosphere sizes larger by a factor of ~ 2. None 
of the upper limits on astrosphere sizes are small enough 
to put the known planets outside of the astrosphere; es¬ 
timates range from 2.9 AU for Kepler-444 to 530 AU for 
Kepler-21. While we calculate an upper limit of 530 AU 
for Kepler-21, this may not be an accurat e estimate as 
the s tar is a subgiant (spectral type F5IV: iHowell et ahl 
l2012t) . The relationship between stellar x-ray flux and 


mass loss found by iWood et al.1 (|2005D upon which our 
estimates ultimately rests was derived for main sequence 
stars, and none of the three subgiants in the sample of 
IWood et al.1 (j2005[ ) agree with the relationship found for 
main sequence stars (although none of the three have 
x-ray fluxes within the range for which this relation is 
a good fit to the data). Additionally, as noted above, 
rapidly rotating stars may have higher stellar wind ve¬ 
locities than slowly rotating stars like the Sun, which 
could affect our estimat ed astrosphere sizes for Kepler- 
448 {vsini = 60 km s~^: iBourrier et ^1120151) . 

The remainder of this discussion is conditional upon 
the Kepler targets that we have investigated actually 
lying within the ISM clouds that we have identified; if 
the stars do not lie within these clouds, then their as- 
trospheres will likely be larger than we have estimated. 
We note again that we only calculate upper limits on 
astrosphere size, conditional upon this assumption. We 
find that Kepler-444, 42, 445, and 20 could have astro- 
spheres significantly smaller than the present-day helio¬ 
sphere; Kepler-32 and 88 (the latter if in Gloud V) could 
have astrospheres similar in size to the present-day helio¬ 
sphere, or smaller; and Kepler -88 (if in Gloud VI) could 
have an astrosphere much larger than the present-day he¬ 
liosphere. Among the habitable zone planet hosts in our 
sample, Kepler-186 has an estimated maximum astro- 
sphere size of 59 AU, somewhat smaller than the present- 
day heliosphere, while Kepler-437 and KOI-4427 could 
host astrospheres much larger than the heliosphere. 

Another system of note is Kepler-444, the first Ke- 
pler planet host that is a member of the thick disk 
(IGamoante eiTalllM^ . As a member of the thick disk, 
Kepler-444 has a large peculiar velocity of 154 km s“^ 
and an old age of 11.2 Gyr. Due to this large velocity 
and old age, we calculate maximum astrosphere sizes of 
2.9 and 4.5 AU if it should lie within Gloud II or III, re¬ 
spectively. This is still well outside the planetary orbits 
in this very compact system; the outermost of the five 
known planets, Kepler-444 f, has an orbital semi-major 
axis of just 0.0811 AU. We note that due to the sys¬ 
tematically large velocities and old ages of the thick disk 
and halo populations, (halo stars can also host planets; 
Kapteyn’s Star may hye a p lanet r fAnsrlada-Escude et al.1 
120141: [Robertson et al.ll201^ . such stars should have sys¬ 
tematically smaller astrospheres than thin disk stars. As 
a result any planets around such stars will be more easily 
descreened by passage through interstellar clouds. Due 
to their high velocities, however, passages through such 
clouds will be shorter than for thin disk stars. 

4.3. Future Prospects 

The distances to most Kepler planet candidate host 
stars are very uncertain or unknown, as these objects 
are too faint to have been observed by Flipparcos. Al¬ 
though Gaia will provide distances to the Kepler hosts, 
until these data are released the ISM clouds in the Ke¬ 
pler field of view can be used as distance markers. Given 
our distance limits on the ISM clouds in the Kepler field 
of view, if absorption at a given velocity is observed in 
the spectrum of a target of interest, this target must be 
at least as distant as the intervening cloud. 

This wo rk will also have bea r ing on a different distance 
indicator. iSilva Aguirre et all (I 2012 D demonstrated that 
they can derive distances from asteroseismology to accu- 
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racies of 5%, but that reddening to the target star is the 
major source of error for more distant targets. Our work 
to map the ISM, the source of reddening, in the Kepler 
held of view will inform the reddening corrections used 
to obtain these distances. 

Gaia will have additional bearing on future work on 
this held, as it will provide more precise distances to 
the Hipparcos stars that we have observed, as well as 
distances to many fainter, more distant early-type stars 
in the Kepler held of view. This will allow us to com¬ 
pile a map of the ISM in the Kepler held of view with 
both higher spatial resolution and extending to greater 
distances. Additionally, it will provide distances to the 
Kepler planet host stars, allowing us to identify a larger 
sample of planet-hosting stars that might lie within ISM 
clouds. 

Further observations are capable of constructing a 
more densely sampled map of the ISM in the Kepler 
held than we presen t here. The Hipparcos re-reduction of 
Ivan LeeuwenI (l2007f) contains more than 300 stars meet¬ 
ing our target selection criteria (but discarding our mag¬ 
nitude limits) that lie in or near the Kepler held of view, 
so there is no lack of potential targets. 

5. SUMMARY AND CONCLUSIONS 

The interstellar environments of planets, especially 
habitable planets, can have an impact upon the plan¬ 
etary habitability and climate. This occurs as the astro- 
sphere size is regulated by the surrounding ISM (the ex- 
oLISM) density and streaming velocity, and in turn regu- 
lates the cosmic ra y hux experienced by the planet (e.g., 
iMiiller et al.ll2006h . Additionally, interstellar gas and 
dust can be deposited directly onto the planet. These ef¬ 
fects can alter the planetary cloud cover, surface temper¬ 
atures, ozo ne layer, surface ultraviolet radiation hux, an d 
more (e.g., lYeghikvan fc Fah'ill2004 IPavloy et ^l2005bl l. 
Thus, as future efforts are made to establish the climate 
and habitability of habitable zone planets, it is important 
to also consider the ambient interstellar environment. 

We have presented the results of a small survey of the 
ISM within the Kepler prime mission field of view, the 
hrst focused on this region of the sky. We have measured 
the Na i and K i absorption towards a sample of early- 
type stars. Using these data we have identified six clouds 
located at distances of less than 450 pc, which we des¬ 
ignate Clouds I through VI. We have found that Cloud 
II must lie at a distance of less than 56 pc, placing it 
firmly inside the Local Bubble. All six clouds likely have 
volume densities greater than that of the LIC. In addi¬ 
tion, we have identified five velocity components which 
are detected along only one line of sight. 

Using the constraints on these ISM clouds, we have 
identified Kepler systems with confirmed planets which 
could lie within these clouds. Using the estimated cloud 
parameters, we have then estimated maximum astro- 
sphere sizes for these systems, conditional upon these 
systems actually lying within these clouds. Most in¬ 
terestingly, we find that the astrosphere surrounding 
the habitable zone planet Kepler-186 f could be smaller 
than that of the Sun {ra^max = 59 AU), while the thick 
disk star and planet host Kepler-444 could have an as¬ 
trosphere just a few AU in size. Additionally, several 
known multiplanet systems (i.e. Kepler-20, Kepler-42, 
and Kepler-445) could have astrospheres much smaller 


than the present-day heliosphere (astropause distances 
of a few tens of AU), while the habitable zone planet 
hosts Kepler-437 and KOT4427 could have astrospheres 
much larger than the present-day heliosphere. We note 
again that these estimates may be significantly in error 
due to uncertainties in the cloud number densities and 
inaccuracies in the astrosphere model. 

While our work has relatively low spatial sampling and 
we can only identify clouds with more than one absorp¬ 
tion component out to a distance of 450 pc, future work 
to observe more early-type stars in the Kepler field of 
view can better constrain the ISM properties in this key 
region of exoplanetary interest. Such a map will also bear 
on a number of other astrophysical applications, such as 
reducing the nncertainties on asteroseismic investigations 
and constraining the distances to Kepler planet hosts. 

Thanks to Katherine Wyman for assistance with the 
ISM fitting routines. Thanks to the staff at McDonald 
Observatory, in particular David Doss, for invaluable as¬ 
sistance during the observing run. M.C.J. is supported 
by a NASA Earth and Space Science Fellowship under 
Grant NNXI2AL59H. 

This paper includes data taken at The McDonald Ob¬ 
servatory of The University of Texas at Austin. This 
research has made use of the SIMBAD database and the 
VizieR catalog access tool, both operated at CDS, Stras¬ 
bourg, France. 
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Table 5 

Confirmed or Validated Planets within 450 pc 


Kepler ID 

KOI 

TeS 

Kp 

d (pc) 

Clouds 


age (Gyr) 

-^planets 

^a,max (AU) 

Ref. 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

(10) 

(11) 


Sun 

5780 

-26.83 

4.8 X 10“® 

Lie 

1.0 

4.6 

8 

154, 121 s 

1 


16 Cyg B 

5772 ± 25 

6.0 

21.21 ± 0.12® 




1 


2 

444 

3158 

5046 ± 74 

8.717 

35.7 ± 1.1® 

II, III 

0.758 ±0.043 

±x.^o_o 99 

5 

2.9, 4.5 

3, 4 

3 

HAT-P-11 

4850 ± 50 

9.2 

36.4 ± 1.3® 

II 


6.5l®° 

1 


5 

42 

961 

3200 ± 65 

15.92 

38.7 ± 6.3'= 

II 

0.13 ±0.05 

> 4.5 

3 

11 

6 

37 

245 

5417 ± 75 

9.701 

~ 66'! 




3 


7 

445 

2704 

3157 ± 60 


~ 90" 

I 

0.18 ±0.04 

> 5 

3 

30 

8 

21 

975 

6131 ± 44 

8.2 

108 ± 10® 

III 

1.340 ±0.010 

2.84 ± 0.34 

1 

530® 

9 

446 

2842 

3359 ± 60 


~ 120" 




3 


8 

410A 

42 

6325 ± 75 

9.364 

132 ± 6.9'^ 




2f 


10 

68 

246 

5793 ± 74 

10.0 

135 ± 10'^ 




3 


11 

438 

3284 

3748 ± 112 

14.467 

145+2°" 




1 


12 

186 

571 

3788 ± 54 

14.6 

151 ± 18" 

I 

0.478 ± 0.055 

2.55" 

5 

59 

13 

10 

72 

5680 ± 91 

10.96 

173 ± 27'^ 




2 


14 

22 

87 

5518 ±44 

11.664 

190 




1 


15 

63 

63 

5576 ± 50 

11.6 

200 ± 15" 




1 


16 

16 

1611 

4450 ± 150'* 

11.762 

~200" 




1 


17 

296 

1422 

3572 ± 80 

15.921 

226+2®" 




5 


12 

1 

TrES-2 

5960 ± 100 

11.3 

~230" 




1 


18 


4427 

3813 ± 112 

15.645 

240l®2c 

V, VI 

u.ozo_o 062 

3.6l?;® 

1 

300, 380 

12 

440 

4087 

4134 ± 154 

15.134 

26lti«" 




1 


12 

441 

4622 

4340 ± 177 

15.142 

284^280 




1 


12 

20 

70 

5455 ± 100 

12.498 

290 ± 30" 

IV, VI 

0.912 ±0.034 

8.81+^ 

5 

54, 63 

19 

32 

952 

3793^®° 

15.913 

303 ± 14" 

I, VI 

0.54 ± 0.02 

3.21" 

6 

180, 99 

20 

2 

HAT-P-7 

6350 ± 80 

10.5 

320ti°" 


1 47+0.08 

2.2 ± 1.0 

1 


21 

421 

1274 

5308 ± 50 

13.354 

320 ± 20" 




1 


22 

45 

254 

3820 ± 90 

15.979 

333 ± 33" 

V, VI 

0.59 ±0.06 

0.69" 

1 


23 

88 

142 

5471 ± 50 

13.113 

339tg" 

V, VI 

0.956l°;°fi 

9 9 + 2.4 

z.z _2 0 

2 f 

140, 270 

24, 25 

442 

4742 

4402 ± 100 

14.976 

342+1°" 




1 


12 

62 

701 

4925 ± 70 

13.75 

368" 




5 


26 

437 

3255 

4551 ± 100 

14.352 

4i7j:^^" 

VI 

0.7071°°®® 

9 q-i-Y .b 

1 

290 

12 

448 

12 

6820 ± 120 

11.353 

426 ± 40" 

V, VI 

1.452 ±0.093 

1.5 ± 0.5 

1 

170, 340 

27 


Note. — Confirmed or validated nearby planetary systems within the Kepler field of view, along with the Sun for comparison. (1) Kepler 
number. (2) KOI number, or alternate name for pve-Kepler planets. (3) Effective temperature of the host star. (4) Host star magnitude in the 
Kepler bandpass. (5) Distance in parsecs. (6) ISM clouds within which the system could lie; see text for more details. (7) Mass of the host star. (8) 
Age of the host star. (9) Number of known planets orbiting the star. (10) Estimated maximum astrosphere size. For stars which overlap with more 
than one cloud, the sizes calculated using the parameters of each cloud are listed in the same order as the clouds in column (6). (11) Reference for 
columns 3, 5, 7, and 8 (unless noted otherwise). A second reference denotes the source of the stellar radial velocity used in computing ra,,max- In 

the interests of br evity, v alues are no t given for columns 7 or 8 if the system cannot lie withi n one of the ISM clouds. _ _ 

References. fl^ |Lang| 1119921') : (2) IWittenmyer et al.l 1120071b (3^ ICa mpante et al.l 1120151b [Nordstrom et ^ II2004I): (5^ IBa kos et all 1120101b 


(6^ IMuirhead et al.l i 20121) : (7) 

Gdli^n^^^^Jj20R b (12) Yh 


IBarclav et al.l l|2013l) : f S) iMuirhead et alj (120151) : (9) j|HOTy^]_et__aL 112 012|b flO) IV^ n Evlen et al.l ||2^ 


(12) [T^ (13) IQuintan^^^aH 


banchis-Oieda et al.|i201^b fl7)IDovle et al.| 1|2011|): (1 8)IO’Donovan et aLHuSnS). , 
II2008I): (22) IKip ping e t al.l 1120141) : f23) IJohnson et al.l 1120121) : (24) INesvorn^ et al.l j201 
IBourrier et al.l il201a). 


(14) Eatalli^^^i 1 1120111) ; 


Gautier "^ "ai 1120121k et al.l 1 I 2 OI 3 I): f21)IPdl et al.l 

: (25) iBarros et al.l 112014) : f26) IBorucki et 1120131) : (27) 


ll2^ 


^ Circumbinary planet. The listed Tgff is that of the primary star; only the mass of the secondary star is known (0.20255 Mq). 
^ Distance from Hipparcos data. 

^ Distance estimated from stellar properties and observed magnitude. 

^ Distance determined from asteroseismic analysis. _ 

® Age estimated using gyrochronology bv lWalkowicz Basril II2013I) . 

^ One transiting planet plus an additional non-transiting companion detected via transit timing v ariations. _ 

s Estimated using our model (first figure; see text). As measured by Voyager 1 (second figure; see lGurnett et £0120131) . 
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